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ABSTRACT 

Immersion lithograidiy is a viable method for continuing the reduction in critical dimension* Much of the improvement 
in image quality in immersion lithography centers around high-NA vector imaging effects and in particular the roles of 
the resist coupling and polarization properties. Electromagnetic scattering from local inhomogenities is considered yet 
emphasis is placed on die unportance of accounting for high-NA vector, immersion and resist standing wave effects. 
A new vector TCC formulation is introduced in SPLAT 6.0. Tlie formulation is tested against a new theoretical 
formula for the Strehl ratio in air at higih-NA and imaging 2D contact patterns at higji NA in the resist stack 
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1. INTRODUCTION 

Immersion lidiography is potentially unpacted by a multitude of effects on many technical fronts. B.J. Lm gave an 
excellent overview of the many factors that potaitially impact immersion lithography in a very early paper [1]. 
Immersion brings benefits such as an option to improve the working resolution, increase the depth of focus or both. 
Techniques for imm^sion lithography to extend the 193 nm node may also be applied to lS7nm using a suitable type of 
liquid [2] to bring continued improvement in linewidth reduction [3]. 

The groundwork for characterizing effects of liquids in optical projection printing is being actively pursued by several 
other groups [2,4-8]. Some of the effects include into^action between the liquid and resist, propagation and scattering in 
inhomogeneous media, high-NA and polarization effects. A number of these effects are also dynamic and can even 
occur on a time-scale within a laser pulse. For example, a projection exposure may produce local pattem dep^d^t 
heating of the resist diat conductively heats the liquid that makes the liquid optically inhomogeneous introducing 
additional lateral scattering. It is important to account for these effects v/hea simulating next generaticm optical 
lithography tools. 

Much of the improvement in image quality in immersicm lithography centers around high-NA vector imaging effects 
and in particular the roles of the resist coupling and polarizaticHi. T. Brunner pointed out the extraordinarily large 
enhancemait of the poor contrast of the TM polarization with the effective high-NA in the resist seen with immersion 
[4]. Smith and Cashmore emphasized the disparity between TE and TM imaging [9]. Various models for high 
numerical q)erture imaging have been proposed. Hagello et. al. first proposed a model for 3D numerical aperture 
imaging applied to optical lithography [10-11] and clearly demonstrated the difference in theTE and TM exposure of 
resist materials. Yeung et. al [12] provide a general vector Hopkins formulation of high-NA effects. They concluded 
that a single term average of the TE and TM effects was suitable up to NA of 0.7 and that resist coupling effects could 
be viewed as an additional aberration term in the pupil. Kirchaer [13] also provides equations for a TE/TM 
decomposition and off-axis rays. Collins [14] deals with an equivalent TCC interaction in three dimensions however it 
is not clear at this point on how the Z^nike polynomial coefficients of much concern in optical lifliography would be 
considered in his model. 

The original intent of this paper was to use simulaticm of both electromagnetic scattering from local inhomogenities and 
high-NA vector imaging to help identify effects in immersion lithography that potentially limit image quality. Section 2 
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briefly surveys the pot^itial for electromagnetic scattCTing from local inhomogenities and explains our decision to focus 
our efforts instead on imaging effects. We quickly move on to a new and fully rigorous reformulation of high-NA 
vector imaging that is suitable for rapid image computation. Section 3 gives an overview of the generalization of 
Hopkins' method of image calculation based on high-NA vector-transmission oross coefficients and the detailed 
derivation is provided in Appendix A. Instead of using a general matrix of nine TCC terms in Eq 21 in Yeung, we 
simplify and consider only six coefficients, i.e. discard the contribution from the z-component as redundant at the 
entrance pupil plane. Several interesting properties of the cross terms and their size versus NA are given in Section 4 as 
well as Appendix B. The transformation conditions under which existing simulators can be utilized to emulate 
immersion lithography systems are given in Section 5. The initial implemoitation of the new hi^-NA vector 
formulation in tiie SPLAT aerial image simulator as version 6.0 is discussed in Section 6. Testing of the new simulator 
is considered in Section 7 and includes a new theoretical formula for die Strehl ratio in air at high-NA. Imaging of lines 
and most importantiy the 2D patterns of contacts are considered in Section 8. 



2. SCATTERING FROM LOCAL IMHOMOGENITIES 

Ibe initial results for inmiersion litiiograirfiy resist images rq)orted by Switkes et. al [2] showed considerable line-edge 
roughness compared to resist images exposed in air. This led the authors to speculate on the possibility of various 
sources of optical noise in the exposure system. Possible scenarios included ghost images from multiple reflections 
back and forth across the liquid, speckle due to scattering from the prism surfaces, scattering from bubbles and 
turbulence and even local heating of the liquid or resist surface during the exposure process. 

The likelihood of these various scenarios were tested or critiqued several ways. M. Switkes provided us a SEM GIF file 
which we analyzed by Fourier spectral methods. No evidence of a ^ost reflection was observed. Further the spectral 
content did not show any evidence of wave contributions from speckle of off-axis waves. The possibility of scattering 
from local heating does not appear to be higfi. This is due to the fact that the heat rise per pulse is at most a few dpgrees 
and dn/dT is small and die fact that thermal difrusion distance in the time between pulses is mudi larger dian a feature. 
Micron sized bubbles were not viewed as a likely source as the hquid surface t^sion would make their internal 
pressure exceed the vapor pressure and the gas would be reabsorbed quickly. 

From die above considerations we concluded that the line-edge roughness resulted from the interaction of the 
immersion liquid widi the resist rather than an electromagnetic effect Meanwhile from the literature and discussions at 
workshops it became increasing clear to us that the primary problem in adueving good image quality in imm^sion or 
any high-NA system is the loss of co-linearity in die IM or radial polarization. We thus focused our energy in 
developing immersion 2-D image simulator including a polarized for eliminating the low contrast radial electric field. 
Since diis simulator includes lens aberrations its development will also allow us at a later date to investigate any macro 
imhomogenities in die immersicm liquid that might be shown exist. 



3. VECTOR IMAGING MODEL 

Our model is, based on gaieralizing die scalar TCCs as first noted by Hopkins [15]. Tliese TCCs describe die 
interaction of light traveling dirough different paths in the optical system before interfering cm the imaging plane. The 
degree of interference is based on the partial coherence sigma of the source. The proposed model builds on the thin 
film formulation by Yeimg [12] akeady in SPLAT as a vector extension to TCCs and follows a similar development as 
Adam [16]. The vector algebra up to the wafer plane results in the same expressions as Flagello [11] when considaing 
die direction of the plane wave strictiy dictated by it's spatial harmonic dius having a predetermined location where die 
wave hits the entrance pupil in our model. Our dgebraic model assumes that Ex and Ey are known on a matfaanatical 
plane just below the mask. These two field components are sufficient to calculate the both the magnitudes and phases 
of a set of plane waves in space that are equivalent anywhere below the mask to the effect of the mask itself. For 
simplicity we assume that these scattered field components do not change with angle of incidence odier than a linear 
phase factor. However, this could be generalized to cover die full vector effects at entrance pupil if necessary by 
considering various zones of illumination over which die source components are invariant. 



We compose our model by tracing one plane wave from the source plane to the image plane in the resist. An overview 
of the formulation will be provided here while the detailed mathematical derivations are provided in Appendix A. 
Figure 1 describes the optical system for projection printing. The column on the left denotes the reference planes in the 
imaging system along with the subscripts used to denote variables plaining to that plane. The column along the right 
lists factors influencing aerial image formation residing at each reference plane. The basic behavior of the optical 
projection printing system is as follows. Light is emitted from the source and travels in a near-normal incident path to 
the mask. The Ught passes through portions of the mask producing a particular near field intensity. The near field can 
be expanded as a smes of plane waves mutual orthogonal in the spatial frequency domain. These plane waves are 
captured at the entrance pupil where the image is reduced by the projection optics. Assuming a diffraction-limited 
system, some of these plane waves cannot propagate out towards the wafer plane. The final aerial image can be 
calculated at any depth in the resist using thin fihn techniques and projecting the E-field mto Cartesian coordinates. 

We now examine the padi of one ray from the source plane to the resist image plane through other planes in the optical 
lithographic system. Without loss of generality, consider ttie path a wave takes whose E-field is oriented in the x- 
direction. Alter the wave scatters from the mask, the electric field hits the entrance pupil and can be split into radial 
and tangential components. Ihe tangoitial component propagates to the exit pupil without a change in orientation 
however die radial oomponoit is tilted partially in the z-direction. Ihe wave reaches the wafer ^ere it is partially 
transmitted into the photoresist depending on the standing wave coefficient from the surface to the image plane. The 
ray changes angle according to Snell's law before readiing the image plane. Accounting for subscripts as listed in 
Figure 1 the representation of the E-field at the image plane is given as: 

ExR= E^t, (0J cosOj, cosip^ - E^t^ S^iO^) sin 0^ (3.1) 
EyR= E^t, Sj^ {OJ sin (p.-^-E^tj^ S^^iOJ cosip^ (3.2) 
EzR=-Vi.^m(^w)sin^;, (3.3) 

where cos0 , U f) _ (g + g^) ^^ ^^^^ 



sin^^=|^*^l e^^cos''[f^) (3.6-8) 



The coordinates (f*, g*) correspond to the pupil position for either the first or second pupil while (f,g) represents the 
point of integration as shown in Figure 2. The obliquity factor, scattering coefficient, and aberrations depend on the 
mask and pupil positions and act equally on all Cartesian componoits. 

In order to utiUze Hopkins method, the TCC must be generalized to the interaction of vector waves oriented at different 
angles \^en passing through the lenis. Since the E-field at the image plane is now a vector quantity, the degree of 
mutual interference corresponds to a complex dot-product operation. Therefore, a vector formulation for the 
transmission aoss coefficient is given by: 

TCc(/',g';/",r)=J„Ji(/',^')(£/-., *4".r)^ (f + f\8 + gj^'(f + f">8 + 8")dfd8 (3.9) 



Figure 2 depicts the region where numerical mtegration of the transmission cross-coefticient J repres^ts the mutual 

coherence circle of radius a, P is the aberration due affecting the first plane wave P*is the aberration affecting the 
second plane wave. The region Q, corresponds to the intersection of three circles: two shifted pupil circles of unit radius 
centered at (f , g') and (f \g") as well as a circle in standard position with radius a. 



4. VECTOR INTERACTION PROPERTIES 

Consider the properties of the vector dot product term added to the transmission cross coefficient. For simplicity of the 
formulation, consider an ideal immersion lens where the full benefit of the refractive index can be utilized. Such a 

situation arises when the rays enter the immersion medium at normal incidence [17]. In this case, ^„ = tj^ for every pair 

of waves. To simpUfy the notation, denote E^^^^, and E^. ^. as Ei and Ej respectively and hide the angular dependence 
of the TE and TM standing wave coefficients in the thin film stack for q=i.,j so that: 

The dot product can be grouped into high NA and low NA terms: 

D{ij) =(£,. •£;p=Ey,owNA{^^p^^.,<Z>p,. -tj)'^E^^^A(0^.O^.(p^ (4.2) 

where: 

Diowna( Oju , , (ppi - (ppj ) = E^^ E^j [Sjj^i Sjj^ J \[cos0a cose^ cos(<p^t -^pih 

+ ^rpi \.^m I ^TE j 

J[cos^jj,sin((^y,-<^pJ] (4.3) 

and 

^0^(0m^^iU><Ppi-^pj'>-E^iErpj[SmiSmj^^^^ (4-4) 

Ihe first term listed in the low NA grouping shows the interaction of the TM ccnnponents ^ose magnitude in the x-y 
plane is proportional to the cosine of the zenith angle for both rays. Additional terms are present at low NA that 
account for differences between standing waves formed from TE versus TM waves. 

The first tam in the high NA equation accounts for the negative interaction of the z components in rays originating 
from opposite sides of the lens. The second term accounts for die alignment of the TE componaits oriented on the x-y 
plane. The second term is die desired component in high NA imagmg. 

One interesting property is the interaction at the image plane of source fields ori^ted in the x-direction with fields 
oriented in the y-direction. In the simplest case whai i.=j, the algebra reduces to: 

D[u^)=a/2)E^E^smi2(p^) [Sm,'J(l- cos'(^, ) (4-5) 

The final multiplying factor accounts for die difference between the TE and TM transmissicm coefficients assuming a 
semi-infinite medium. Figure 3 displays the effect of this cross term wiien used in a dry system versus in an unmersion 

system relative to the odier coefficients in the dot product, i.e. only (1- COS^ ~" ^, ) ) is compared. We assumed a 

refractive index of 1.72 for the photoresist and a refractive index of water of 1.437 at 193 nm [9]. 

In the dry imaging systems, this term has a magnitude of approximately 5% relative to the self-interaction terms in a 
system with a NA of 0.6S when high NA imaging effects start to appear. At Brewster's angle in the dry system, diis 
cross-term exceeds 10% of the relative magnitude of self-interaction terms as NA approaches 0.85. In comparison, the 
immersion system reaches 5% cross-term interaction only as NA increases beyond NA of 0.92, well beyond the 
Brewster Angle for the water-resist interface corresponding to NA of 0.762. 

Appendix B illustrates additional basic properties of the dot product such as algebraic substitution of x and y oriented 
E-fields, verifying that a self-intaraction of unity occurs in die on-axis case and that rotational symmetry exists in die x- 
y plane. 



5. IMMERSING EXISTING SIMULATORS 



This section outlines a set of four substitutions fliat allows any aerial simulator correctiy account for the effects of 
introducing an immersion liquid between the exit pupil and the wafer plane. Several assumptions before making made 
Miien proposing these changes and they should not be applied blindly without understanding their in:q)lications to the 
particular simulation program in use. The modifications assume that the refractive index of the immersion liquid is 
relatively low ( n < 2 ) so that when the demagnification factor is scaled the rays at the edge of the entrance pupil still 
have a negligible component in the z-direction. Otherwise, a more rigorous model will be needed to model the 3D 
orientation of the x-oriented and y-oriented E-fields at the entrance pupil. Lasdy, it is important to understand the 
implications of these changes in other parts of the program as the following formulation is based on the manipulation of 
internal variables m SPLAT 5.0. 

Introducing an immersion liquid alters the imaging wavelength, the zenith angle of incidence on the wafer and also 
alters the E-field coupling into the resist thin film stack. Assume that NAs is die effective numerical aperture in air that 
the user provided for the aerial image simulator. The first obvious substitution is reducing the wavelength by nmed. In 
order to correct for the numerical aperture, using Snell's law and the definition of numerical aperture: 



To maintain the E-field coupling into the resist thin film stack, consider Snell's law and other related fimctions: 



Tlie reflection and transmission coefficients for parallel and saggital orientations can be expressed as functions of 
elementary sine and cosine relations of the incident and transmitted zenith angles. Therefore, the third scaling rule is to 

reduce the indices of refraction of all materials in die optical padi by a factor of . The fourth scaling factor refers 

to the entrance pupil acceptance cone. If NA represents the effective numerical aperture in air, the numerical aperture 
at the entrance pupil is NA/Mdemag \^ere Mdemag denotes the optical demagnification factor. In order to maintain the 
numerical aperture at die entrance pupil ^en scaling NA to NA dira die optical demagnification factor must be 
reduced by a factor of . This adjusts the obliquity factor in an appropriate manner. 

As an example, consider Figure 4 depicting the entrance and exit pupils of a 4x demagnification system with NAo=0.8 
at the exit pupil for a dry imaging system and an ideal wet imaging system immersed with a liquid of n=l.S to achieve 
an effective NA© of 1.2 in air, Assmning that the optical system is diffraction limited at the exit pupil, the physical ray 
cone dictating the numerical aperture is always constrained to 0.8. However, the zenith angle from the leas to the 
immersion liquid at die exit pupil can be larger since there is less constraint on the critical angle leading to total internal 
reflection. In effect, this allows a larger numerical aperture of rays at the entrance pupil to be transferred through the 
exit pupil. The increase is proportional to the refractive index of the immersion liquid. Therefore die demagnification 

factor should be decreased by a factor of for the purposes of calculating die obliquity factor and normalization 

constant as detailed in Cole et. al. [18]. Image reconstruction in the immersion liquid requires that spatial frequencies 
be scaled by Umed equivalent to a reduction in wavelength. In summary, the following transformations can be used in 
order to simulate immersion lithography: 



sin fff _ NAf _ 
sin 0i NA^ ftf 




(5.1) 



sin^j = — sin^. 
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6. IMPLEMENTATION 



Using Hopkins theory of partially coherently unaging, the total image is the integration of the transmission cross- 
coefficients over all available pairs of object spectra 0(f,g) v^ere the scalar TCC has been replaced with Eq (3.9) as 
described in Section 4 and derived m Appendix A: 

1= jj llTCC(f',8'-,r,g'')d (/',g')d'{/",8") 

.X tx^[-iiK[{r-nxHg'-r)y)4'dr drdr (6.i) 

Equation (6.1) is the basis for SPLAT version 6.0, which accounts for the vector nature of the interaction of light The 
current version of this code has been implemented in Fortran within die existing program ardiitecture of Splat 5.0. The 
program is currently in a non-optimal state and Ongoing verification of this image model is in progress by comparing 
SPLAT 6.0 to existing results in works such as [4],[7],[9], [19]. Algorithms for the detection of symmetry were 
bypassed to ensure correctness leading to an increase of up to 16 times as many TCC calculations involving numaical 
integration. In the case of calculating the image in air, no additional loops are required so only the symmetry penalty is 
incurred. In die case ^en the image plane is inside the resist, there has not been an appropriate look-up table for the 
vector case so there is an extreme penalty for resist-coupled cases at the moment This will likely improve in the near 
future as functionality is added to SPLAT 6.0 and the symmetry detection algorithm is generalized for the vector TCC. 
In the optimal case, die incurred p^alty could readi a factor of 2 or 3 over die scalar case at hig}i NA even for image 
plane calculations inside the resist 



7. VERIFYING THE VECTOR MODEL 

This section formulates and compares die calculated and simulated values for a small 2D feature in SPLAT 6.0. We 
position an isolated 0.2 UNA contact at the center of the mask. This feature is sub-printable and approximates a 
pinhole so we assume that die mask scattm light in all directions, filling die entrance pupil uniformly. In order to 
simplify the TCC integration, we tested our vector model with a partial coher^ice of O.OS. Hie transmission coefficient 
between die exit pupil and immersion liquid is identical for all pupil positions assuming an ideal immersion system and 
can be factored from the expressicms in Eq (3.9). Hierefore, die E-field in air vAiea neglecting the effects of the resist 
is: 

ExR= - sin (p^ (7.1) 

EyR= E^ cos0^ sin (p^ + E^ cos(^^ (7.2) 
EzR=-£^sin^^ (7.3) 

Supposing now that the E-field incidait on die mask is oriented in die x-direction, we can substitute the radial and 
tangential fields. Ilie E-fields can be integrated over the pupU plane fi*om 0<^p <2n and 0 < r < 1. Integration of EyR 
and EzR over the pupil plane yield zero but over Exr: 

ExR= E^[ J jcosff^cos^ (p^rdrde- J j sia^ (p^rdrdO] (7.4) 



Noticing that die zenith angle is a function of r and is constrained fi:om 0<6j^< NA. Using sin(6R)=rNA and 
substituting p=rNA, integration yields: 



2(l-(l-iVA^)^^^) 

JJ ^ ^ ^ 

This can be viewed as a perturbation of the scalar case as NA changes so that the normalized field is: 

£^=1-A where A= (7.6) 

A comparison of the scalar versus vector versions of Splat against equation (7.6) is shown in Figure S. in the 

scalar versim 5.0 remains constant ^ereas the vector version 6.0 predicts the value to within the numerical oror of the 
simulation program. Be aware that this case does not validate diat the vector model is correct, it has just shown that 
reasonable assumptions and careful programming were used to verify a simple case. The derivation was carried out for 
a demagnification factor of one where there are non-negligible z-components. We are developing additional scenarios 
>^ere analytical analysis results in a closed form solution that can be used to test our vector model. 
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8. SIMPLE PATTERNS 

This section shows three simple cases ^ere the new models in SPLAT 6.0 can be used. The first case shows: the 
difference between scalar and vector imaging at hi^ NA. The second case compares image formation within the resist 
for a dry and immersion system with comparable parameters. The third case displays the effect of introducing a pupil 
plane filter and imaging solely with die TB or TM waves. 

We chose to use a 0.55 X/NA L=S PSM patten to compare the use of scalar and vector imaging predictions at higji NA. 
The optical system we used was 4x demagnification, sigma =0.1 and NA ranging from 0.7 to 0,9. Photoresist 
interaction was neglected for these simulations. Figure 6 and Figure 7 plot the intensity in air for the Splat 5.0 and 
Splat 6.0 respectively. Both simulators predict the fi-equency-doubling ph^omenon however Splat 5.0 overestimates 
the peak intensity that should theoretically be 2.0. Splat 6.0 predicts that flie quality of the nulls will eventually be 
degraded since there is no longer total vector cancellation at hig^ angles of incidence. 

For the second case of simple patterns we considered an isolated contact of size 0.7 UNA at NA=0.85 and sigma of 0.3. 
A homogeneous photoresist film 30Qnm tiiick with n=1.72 is used choosing a substrate producing a 15% E-field 
reflection coefficient. The aerial image is calculated inside the resist at a depth of 50 nm. The aerial image in the resist 
using 157 nm in air is shown in Figure 8 while the same intiage in the resist using using 193 nm in water immersion is 
shown in Hgure 9. Both die scalar and vector models show a greater degree of coupling to the resist stack vAim using 
water immersion as a result of a higher transmission coefficient into the resist The scalar model overestimates the 
image formed using the vector model in Figure 8. The vector model shows that the resist coupling at 193 nm as 
mitigated by immersion gives nearly a 4x improvement in the brightness of smaU 2D features. 

Finally we observe the predicted effects of inserting a polarization filter at the pupil plane [20] to improve image 
quality at high NA. We use the same mask topology, imaging parameters and resist image plane as above except we 
use 193 nm in air. The intensity cut lines for the TE image and TM image at a 50 nm depth inside the resist are plotted 
in Figure 10. The TE and TM intensities can be directiy compared to those shown in Figures 8 and 9.The image formed 
by filtering TM waves at the pupil plane has higher peak intensity and better image quality than its counterpart The 
major disadvantage in impl^enting such a filter is diat approximately 40% of the non-filtered intaisity is imaged in 
the resist raising issues such as insufficient exposure, thermal absorpticm or reflection by a real filter. 



9. CONCLUSIONS 

We introduced a new high-NA vector transmission cross-coefficient The formulation permitted us to separate the 
high'NA and low-NA terms, present qualitative insight and show that orthogonal source x and y e-fields lose thdr 
orthogonality more severely in dry systems compared to immersion systems. We also introduced and verified a new 



fteoretical formula for tbe Strehl ratio in air at high-NA using SPLAT 6.0. We considered imaging of 2D contact 
patterns in Section 8. The peak intensity for these small features in the resist improved fourfold with inmiersion 
primarily from improving the resist coupling. 
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APPENDIX A 



This section provides the method used to derive the transfer of E-fields from the source plane to the image plane inside 
the resist Figure 1 is used as a reference for defining the variables and subscripts for our vector model. Assume that 
each linearly polarized plane waves hit the mask whose polarizations are either in the x-direction or y-direction. Each 
of these waves has an electric field that can be split into a radial component with subscript r and a saggital component 
with subscript <!> depending on where the wave hits the pupil: 

x-oriented: = £^ COS (p^ = -E^ sin ip^ (A.1-2) 

y-oriented: E^ = Ey sin E^ = Ey cos (p^ (A.3-4) 

Analysis will proceed to the exit pupil where all rays transition between the exit lens and the medium above the wafer. 
Aberrations based on the OPD value on the pupil map affect all fields equally contributing to a temporal phase shift 
When the rays enter the medium at the exit pupil, their amplitudes are modified by the parallel and perpendicular 
transmission coefficients. Hie lens designer would supply infonnation about the ^trance angle versus pupil position. 
To simplify the problem, assume that the rays enter Ae milium at nearly normal incidence for all pupil positions across 

the lens. This is the ideal inmiersion case where the imaging resolution increases by . The transmission 

coefficient for either ray reduces to: 

2n, 

h=h = so that E^ = E^t^ and E^ = E^t^ (A.5-6) 

where nias is the refractive index of the lens. This expression is constant for all rays and can be factored outside of the 
TCC integration. Combining these coefficients with the standing wave coefficients is possible depending on the type of 
prism used to couple light into the imma-sion liquid The tangratial field emerges from the exit pupil without any axial 
tilt however the radially oriented field is tilted at an angle Ow with respect to die z-axis depending on its location on the 
pupil: 

COSd„=.f^ ,k^=^,k^=^,k^ =^k]+kl (A.7-10) 

Once the wave hits die resist, it goes through a multitude of reflections and transmissions modeled by upward and 
downward propagating plane waves derived by Yeung [12]. These coefficients are different for TE than for TM waves, 
corresponding to ^e tang^tial and radially oriented fields used in the derivation tiius far. Denote these standing wave 
coefficients by Stc(9w) and Stm(Bw) so that: 

ErK=SmWE^ and 8^,(0 JE^ (A.11-12) 
The following angular relations hold based on die orientation of the problem space and Snell's law: 

sin 9^ = sin ^j, = = <l>p (A.13-14) 

Finally, the field is decomposed into its x, y and z coordinate orioitations so diat: 

^= ^rR cos^j cos^j - sin <f>g^ (A.15) 

EyR= £rff cos^jj sin^>jj+ E^ cos^^ (A. 16) 

EzR=-£rfjsin^jj (A.17) 

Substituting the results from Equations (AS) through (A.14) into (A.15-17) yields: 

ExR= E^t^ (OJ cose^ cos<pp - E^t^ Sj^(9J sin (p^ (A18) 

EyR= E^t, 5™ (OJ cosO^ sin + E^t^ S^idJ cosip^ (A19) 

EzR= - E^t^ Sjj, {e„ ) sin (A.20) 



where cos^„.-7= ^f'^f^ dnM (8 + 8 ) (A.21-22) 

' ^|(f+fy+(8+8y yiif+rr+is+gy 



COS0g= Jl- 



^med J. 



e„ = cos"* (yjl-k^) (A.23-25) 



The point (f.g*) used in the calculation of the azimudial angles is the pupil position for either unit circle in Figure 2 
centered at (f,g') or (f'.g^'). 

APFENDKB 

This section elaborates on Appendix A by simplifying the dot product expression and considering some vector 
interaction properties under different circumstances. To simpUfy the notation, doiote E^. ^. and Ey. as Ej and Ej 
respectivdy. Fen- two rays i. and j, the dot product expression using Eqs. (A.18-20) becomes: 

£j = ExRi ExRj + EyRi EyRj + EzRiEzRj (B.l) 

= l^nt^K Sm(^m) cosO^ cos<ppi- E^.t^ sin^zJ^,.] 
^ ^^>pAi ^mi^y^) cosOgj cos(ppj- E^jhj Smi^wj) sin^^py] 
+ ^E^ASm(^wi) cosej^sm<l)pi+ E^it^SjEiOyn) cos^^p,.] 
(O^j) cos0gf sin^p,.+ E^jtj^j Sj^iO^j) cos^^,] 
+ lE^f«Sj„(.^M) sm0„]^[E^jhjSj„(0^) sm0^] 

where x is used for scalar multiplication. For simplicity, factor out (he transmission coefficients and simplify the 
notation for q = i. ot j by denoting: 

Group terms in Eq (B.l) by field component products at die pupil plane reduces die dot product to: 

(E,»Ej)= {B.5) 

The dot product is a function of the difference between pupil angles as expected and exhibits rotational symmetry 
across the pupil. First, consider X wave or Y wave self-interaction. In the simple scenario when standing wave effects 
are neglected and assuming that i. and j are equal then tiie dot product of two E-field vectors in air should yield the 
original field value squared. From (B,5): 

E,*E,=El, + £^,£^,(0-0)+ V= Kt + V <B-6) 
Fot x-waves and for y-waves, use (B.5) to get: 

El, + V= ^/cosVp,+^x'smVpi=^/ (B.7) 

El, +£^/= sin^ <P^,+e; cos' <p^,=E^' (B.8) 



Second, consider X and Y wave interaction when i = j. y/hea wave i is transformed into an x-wave and wave j is 
transformed into a y-wave, the interaction whai i=j should be zero. Use (B.5), drop E product terms then substitute x 
and y waves from (A. 1-2) to get 

{Ei9Ej)=E^ Eyj cosippi siniffpj [Sna Sj^jl-E^ E^ sin^^,. cos^^^. [Sj^i Sj^j\ (B.9) 

When i,=j then, 



( 



E,^Ej)= i E^ E^ sin(2«>^,)( b^,' J ■ [Sj^'\) 



(B.10) 



The conclusion is that x-waves and y- waves at the same point in the pupil plane don't show orthogonality at the wafer 
plane as NA increases. In order to gauge when these interactions are significant, consider a semi-infinite medium 
below the wafer plane. The ratio of the transmission coefficients for the two polarizations is: 



2sin^,cos^, 



2sindf cosd^ 



so that. 



b™/J- [57^/])= bm.'Jd- cos^(^, 



(B.11) 



(B.12) 
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Figuie 1 Optical Factors affecting a Projectioii Printing System Figure 2 Region of Interaction Leading to the TCC 
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Figure 4 Giange in Entrance and Exit PupO Cones when adding 
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Figure 5 Comparison of Theoretical vs Computed Aerial 
Imag& formation of a 0.2 ^NA pinhole in SPLAT 5.0 and 
SPLAT 6.0 
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Figure 6 Aerial Image in Air, Scalar Model, 0.55 }JNA L=S PSM 
with increasing NA 
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Figure 7 Aerial Image in Air, Vector Model, 0.55 l/NA L=S Figure 8 Air and 157 nm tool aerial imag in resist @ z=50 nm, 
PSM with increasing NA NA^.85, 0.7 X/NA isolated contact 
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Figure 9 Water and 193 nm tool aoial Image using in resist @ Figure 10 Aerial images through pupil plane filter, measured in 
2=;50 nm, NA=0.85, 0.7 7JNA isolated contact resist @ z?=50nni, NA=0.85, 0.7 UNA isolated contact 



